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A3STRACT 

OPTIMIZATION OF HYBRID EVAPORATIVE COOLING AND 
AIRCONDITIONING SYSTEMS - AN ECONOMIC APPROACH 


A Thesis Submitted 
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for the Degree of 


MASTER OF TECHNOLOGY 

by 

K. SEKAR 
to the 


Department of Mechanical Engineering 
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May , 1986 


A generalised computer programme has been developed 
for the cooling load calculations on hourly basis and the 
sol-air temperature variation for different orientation of 
walls. The effect of heat capacity and timelag factor of 
the structures on the cooling load estimation has also 
been studied. 

The single-stage vapour-compression airconditioning 
system has been optimised to get the minimum total cost, 
for refrigerants R— 12 and R-22. This is compared with a 


( xiii) 


hybrid system of evaporative cooling and airconditioning as 
to what extent, the latter is more economical than the 
former. The optimum results have been found out for the 
inside design conditions of 24°C, 65%; 27.5°C, 56% and 30°C, 
60% (dry bulb temperature and relative humidity) beirsg 
obtained on the basis of Fanger's comfort equation £6] . 
Various combinations of the operating parameters are also 


studied 
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CHAPTER 1 
INTRODUCTION 


1.1 A IRCOND IT ION INC SYSTEM; 

Airconditioning is constantly absorbing and apply- 
ing the latest discoveries of science and engineering 
in an effort to ease man's environmental problems on 
earth and in space. The term, airconditioning means the 
act of putting the air in the proper condition for the 
desired use. This may mean cleanliness, the temperature 
and humidity of the air. Airconditioning for human 
comfort is becoming almost a prerequisite for success in 
the merchandising field in all except the cooler regions 
of the country* The gain in customer’s good-will is 
reflected by such an increase in business that the air- 
conditioning pays more than for itself. As any other 
engineering system, the airconditioning system is controlled 
by various system parameters. To achieve maximum utilisation 
and effective implementation, an air conditioning system 
has to be economical. In view of this, the operating 
parameters of the system are optimised on the basis of 
cost as well as functional energy constraints whichever 
is contemplated. To carry out this objective, there are 
alternatives available, namely modification of an existing 
system and design of a new system. The present work is 
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done on the basis of the formes* 

1.2- SELECTION OF ENVIRONMENTAL ZONES AND OUTDOOR 

TEMPERATURES s 

Eventhough, there are many app locations , the 
primary function of an airconditioner unit is considered 
to create temperatures artificially which will be comfor- 
table for human beings. Many investigators have done 
extensive work on temperature zones for human comfort. 

Fanger [ 6 ] elaborated the methods of getting comfort 
conditions and also suggested a comfort equation after 
carrying out extensive field studies. In India, Ramamoorthi 
[15] has analysed the equation for the Indian weather 
conditions after taking into account of several variables 
and suggested the comfort conditions, namely inside dry- 
bulb temperature, relative humidity and velocity of air* 

Singh and Prakash [14 ] presented a new ASHRAE comfort chart and 
recommended a region of design conditions for the comfort 
environment. The result of the field studies carried 
out by Malhotra [10 '] for the conducive atmosphere pertain- 
ing to Indian environment has also been taken into consi- 
deration. Based on the energy cost, Whitner [18] has got 
the higher comfort conditions. Several investigators [1 ] 
assumed the daily to outdoor temperature variation as a 
sinusoidal function which results in Fourier series. In 
India, Ramamoorthi [15] suggested that the air conditioning 


system can be substituted by evaporative cooling system, 
whenever feasible* 

The design of an aircond itioning plant is mainly 
dependent on the cooling load , a function of indoor and 
outdoor temperatures and internal energy release. The 
indoor temperatures are chosen on the basis of comfort 
industrial requirement. The size of the building and 
various capacities of the equipments are also the primary 
factors influencing the cooling load calculations. So, 
great care should be taken in proper selection of the 
equipments to get the total cost as minimum. It has been 
a common practice to use air-cooled condensing unit for 
both residential and commercial applications. This 
dispenses with water piping, cooling water tower and 
water treatment, thereby eliminating their costs. But the 
head pressure »f the air-cooled system is usually 10% to 
20% higher than that of the water cooled system. Hence 
the air-cooled system is not used in case of large cooling 
requirement. For smaller units /the cost and conveniences 
dominate over other factors. Hence, the air-cooled 
system is considered in the present analysis. 

1 . 2.1 ANALYSIS OF COOLING LOAD CALCULATIONS } 

There are number of parameters affecting cooling 

load calculations, often difficult to define precisely, and 

i 

always intricately interrelated. Most of these parameters 
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can be assembled under the classification of space heat 
gain. This is necessary because different fundamental 
principles and equations are used to calculate different 
modes of energy transfer. The major parts of the space 
heat gain are (a) heat transfer through the structure of 
the building/ (b) heat generated within space by occupants/ 
lights and appliances and (c) ventilation and infiltration 
of outdoor air • For the evaluation of structural heat 
gain/ there are different techniques available/ namely 
(a) sol— air temperature method, (b) cooling load temper- 
ature difference method, (c) finite difference method and 
(d) finite element method. Among these, sol-air temper- 
ature technique is mostly preferred due to its simple 
algorithm. The cooling load temperature difference 
method [ 3] is implemented by using the transfer function 
method to compute the one— dimensional transient heat 
flow through the building structure. Kadambi and 
Hutchinson [9 ]have described the finite difference method 
applied for the heat transfer through the walls. Huebner [7] 
has elaborated the finite element method in the applica- 
tion of heat flow through structures. The thermal capacity 
and thickness of the structure are interrelated and 
presented in the graphical form [ 2 ] . The relation between 
the infiltration load and the volume of the conditioned 
space is explained in [3,4,12 ] « Woods et al [19] have 
modified the ventilation standards by giving emphasis on 


energy saving . The heat generation within the conditioned 
space has been elaborately given in [ 12 ] talcing into 
account of metabolic ratings. The procedure for hourly 
cooling load calculation is elaborated in [3 ] . 

1.2.2 METHOD OF OPTIMISATION: 

In genera 1, optimisation can be defined as the act 
of obtaining the best result under the given circumstances. 
In design/ construction and maintenance of any engineering 
system, engineers have to take many technological decisions. 
The ultimate goal of all such decisions is to either 
minimize the effort required or maximise the desired 
benefit . Since the required effort or desirable benefit in 
any practical situation can be expressed as a function of 
certain design variables, optimisation, in particular, 
may be considered as the process of finding the conditions 
that results in the maximum or minimum value of the 
function, known as objective function. It is very rare 
that a practical design problem to be unconstrained. 
Considering this fact, if the expression for the objective 
function and the constraints are fairly simple in terms 
of design variables, the classical methods of optimisation 
technique can be used to solve the problem. On the other 
hand, if the problem involves the objective function and 
the constraints which are not stated as explicit functions 
of design variables or which are too complicated to 
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manipulate, we should go in for some specialised numerical 
methods of optimisation. The basic philosophy of these 
methods to produce a sequence of improved approximations 
to the optimum value. In the present work, one of such 
methods called as 'Penalty Function Method' [16] is 
adopted to obtain the optimum design variables of the 
airconditioning system for the minimum total cost. The 
application of this iterative procedure is made possible 
due to the availability of the computer. These methods 
are clearly elaborated in [16] . 

1.3 PRESENT STUDY; 

In the present study, a generalised computer 
programme has been developed for the accurate evaluation 
of the cooling load considering daily outdoor temperature 
variation based on actual as well as approximate. In 
this regard, the effect of timelag factor and heat 
capacity of the structures on the structural load has also 
been found out. The cooling loads have been calculated 
on the basis of actual hourly daily temperature variations 
and the sinusoidal temperature variation based on maximum 
and minimum temperatures. The multiplier has been developed 
to facilitate the utilisation of the suggested sinusoidal 
temperature variation; because these data are available 
for most of Indian cities [11]* 
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For the iterative optimisation technique, a 
versatile computer programme has been developed which 
takes care of various factors involved in the total cost 
of a vapour compression system. Even though there are 
large number of operating variables present in an air 
conditioning system, the most important factors have been 
selected as design variables, giving emphasis on the 
feasibility of control. For these variables, the constraints 
are formulated on the basis of practical consideration* 

A comparative study has also been carried out between the 
airconditioning systems operating with and without evaporative 
cooling, for the duration, April to June, in terms of total 
cost. This aspect is studied for both R-12 and R-22 
refrigerants. 
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CHAPTER 2 

COOLING LOAD ESTIMATION 


2.1 INTRODUCTION: 

The design of an airconditioning plant for a 
building depends on the use of the building. The cooling 
load calculation needs the size of the building/ its 
structural details and orientation. In addition/ the 
following quantities are considered: 

i) Design conditions: Indoor and outdoor conditions/ 
ii) Instantaneous heat load: Sensible and latent. 

For the given inside design condition/ the cooling loads 
are calculated as follows: 

Sensible cooling load: It comprises the heat 
transfer through structure including solar heat gain. 

The sensible heat part of the infiltration and ventilation 
is calculated from the known amount of the air and subsequent 
difference between enthalpies of the inside and outside 
air. The sensible heat part of the heat release from 
the occupants is calculated from the standard data 
available giving heat release from human beings/ ratings 
of fans# lights/ etc. 

The latent heat part of the cooling load hag been 
calculated from the infiltration and ventilation air/ 
latent heat release from the occupants inside the room. 
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2.2 SELECTION OP INDOOR CONDITIONS: 

The variation in the weather, building occupancy, 
and other factors, affecting the load, necessitate 
carefully coordinated controls to regulate simultaneously 
the components and equipments to maintain the desired inside 
conditions. The design indoor conditions play a key role 
in the design of any comfort airconditioning system. The 
field studies [lo], reveal that in India, for hot humid 
climates, the comfort zone is between 22°C and 24.5°C, 
effective temperatures. The analysis [15] carried out on 
the basis of Fanger’s comfort equation recommends the 
comfort conditions for Indian people as T^ b = 24 °C, 

RH = 0.65, V = 0.8 m/min; T db = 27.5°c, RH = 0,56, 

V = 1.1 m/s, and T^ b «= '30°C, RH = 0.60, V = 0.7 m/s. 
Since, these results lie in the comfort zones recommended 
in [5, 6,14], in the present analysis the same combinations 
are taken into consideration. 

2.3 HOURLY COOLING LOAD: 

Most of the factors of the cooling load vary in 
magnitude over a wide range during a twenty four hour 
period. As the cyclic changes in the load components are 
not usually in phase with each other, the analysis is 
required to establish the resultant maximum cooling load 
for a building. 
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2.4 HOURLY VARIATION OF OUTDOOR CONDITIONS: 

Hourly variation in the outdoor conditions is 
available for the limited cities. But the cooling load 
depends greatly on the ambient temperature variations 
as compared to that obtained on the basis of design 
conditions. So a mathematic model [1] has been adopted 
which gives a close approximation to an average value 
equal to the mean daily temperature plus an oscillating 
component with an amplitude equal to one half the mean 
daily range. For most practical applications/ this can 
further be approximated by a constant with a superimposed 
sine wave having a 24 -hour period. So, the final form 
of the- ■expro&siQjQ. comes out to be a harmonic series of 
first order: 


T = A + B Cos ( 15 t - C) 
o 


( 2 . 1 ) 


where T 


= outside air temperature at any time, °C 
o 

t = time in hours 

A,B,C ss the constants being obtained by applying 
the boundary conditions. 


dT 


dt 


= 0 when 


f 1 


max 


( 2 . 2 ) 


it = t . 

T . occurs at (Sunrise-1) °hrs? tJ*" occurs 12 hrs thereafter, 
mm ra ®-* 


So, the above model is solved for the autdoor air temperature 
variation by knowing the minimum and maximum temperatures 
with their corresponding timings, in a day. The result of 
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this equation has also been verified with actual variation 
of the temperature and the comparative study has been 
done in a graphical form. Fig. 4.1, for Kanpur. For this 
purpose, the data for the actual variation of the temper- 
ature and relative humidity for the period 1982, 1983 
and 1985 have been procured [8] and given in tabular form 
in Appendix -A. 

2.5 SOLAR RADIATION: 

The solar radiation is the sole factor for the 
summer air conditioning system. The total incident solar 
radiation is given by: 

I t = I pN CosQ + I d , k W/m 2 (2.3) 

where the product of 1^, direct normal radiation and the 
cosine of the incident angle 0, stands for the irradiation 
on a surface from the sun and 1^ refers to the diffused 
radiation. Simplified general relations for I DN and 
for a clear sky are given approximately by, 

2 

Idn = A i/ e;x P (Bj/SlnP), k W/m 
- C 1 ‘ k H/ “ 2 

where, 

A-j_ = apparant solar irradiation at an air 

2 

mass = 0, k W/m 


2.4( a] 
2.4(b) 
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= atmospheric extinction coefficient , 

= diffuse radiation factor, 

F = the angle factor between the surface and 

o o 

the sky, 

= 0.5 for vertical surfaces, 

= 1.0 for horizontal surfaces, 

= ( 1 + Cos0r_)/2 for any other inclined surfaces 
and (3 = the altitude angle. 

The values of A^ , and of [3 ] are given in the 
Appendix -B. 

2.5.1 ESTIMATION OF SOL-AIR TEMPERATURE: 

For heat transmission calculations, it is convenient 

to combine the effects of outside air temperature and 

solar radiation clustered into a single fictitious quantity 

called sol-air temperature (T 1 ). The rate of heat 

sol 

transfer, q, from outside to inside surfaces of a sunlit 
structure may be written as, 

q = h (T - T J + a. . I* (2.5) 

^ OO W,0 t O 

where, 

h = outside air convective heat transfer 
o 

2 

coefficient, k W/m - °C 

T = outside air temperature, °C 
o 

T = temperature of the outside surface, °C 
w,o 
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= absorptivity of the structure/ 
and I = total intensity 01 solar radiation/ k W/m * 

The total intensity of radiation is calculated by the 
primary angles namely/ the latitude (1)/ hour angle (a 
and the sun's declination/ d. Besides this/ there are 
several other angles as the sun's zenith angle x , 
altitude angle (3 , and azimuth angle Y . In Appendi:x-C, 
the equations for the estimation of fchese angles are 
given. 

The equation (2.5) can be modified as 

q = h (T . - T ) 

M O SOl W/O 


where , 


sol 


= T 


a t •‘ I c. //h o 


The outside air convection heat transfer coefficient, 
hl^is calculated as follows: 


= 14. 27 8+3. 5 59 .V w 

for 

very smooth surfaces. 

2.7( a) 

= 18. 4 23+3. 81. V w 

i or 

smooth surfaces. 

2.7(b) 

= 26. 55+5. 0663. V w 

for 

rough surfaces. 

2.7(c) 

= 28. 64+6. 364. V w 

for 

very rough surfaces. 

2.7(d) 


where V is the design velocity of the outside air in 
w 

Km/h. In the present study, the variations in the sol- 
air temperature for a day, on different orientation of the 
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surfaces presented in a graphical form,, Fig. 4. 2. 


2.6 HEAT TRANSFER THROUGH STRUCTURES: 

A fraction of the total heat energy is absorbed 
by the structural material and there is also a timelag 
for the heat transfer from the outside surface to the 
inside surface of the structure . The thermal capacity 
of the walls and ceilings, c tll ,is taken as. 


c th = m. c p = P . A . x. c p (2.8) 

where rn, c^, P , A and x are mass, specific heat, density, 
cross-sectional area and thickness of the structure, 
respectively. 

Taking into account of timelag and decrement factor ,ti) 
the actual heat transfer through the structure at any 
time l t' is given by 


•t S [ A j lT om - V + Vr A j (T o(t-x> - T om> 1 


0 , = 


j=l 


( 2.9) 


where T om , A- and "t are the mean sol-air temperature, 
decrement factor and timelag factor, respectively. In the 
above expression, U, A, and t) re ^ er oVeI,a l-^ 

heat transfer coefficient, cross sectional area and 
temperatii&e at time (t— t), respectively. Here summation 
symbol stands for four walls and the ceiling. 



K) = 0.00072 kW/m°C 
K 2 = 0.00130kW/m°C 


Fig. 2.1(a) Structural details of the wall. 


K 3 = 0.00072 kW/m°C 
K 4 = 0.00173kW/m°C 

K3 

k 3 

Fig. 2.1(b) Structural details of the ceiling. 


MORTAR 


do T 0l ho 


f 1** T? ^ \ * ;■> o/-- • VJ 
C ■. 1 • •' V } ♦ *• .'•'** ., • /I *. - \ • < . Ly .* jj d 

: - \ 

> ! 1 11 t • • • 1 • •?.! 


CONCRETE 


^ ^ ^ ^ ^ ^ ^ r ^ ^ ^ ^ ^ ^ 

I l M i qi T i.- h 
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Referring to the Fig. 2.1, the heat conduction 
through structure can be written as, 


( T 


w /°. 


“ T. . , ) 


q . = h (T - T ) = . _ 

st o o w,o 1 =3 


W L 1 


S (x ± /K,) 

i_l 


h . ( T . -T . ) 
l w ,x i 


( 2.10 


The above equation can also be written as. 


st 


= U 


(T, 


- T i> 


( 2.11 


with 


U 



i=3 

S 

i=l 


>: . 
i 

K. 

x 



T q = outside air temperature, °C 
T^ = indoor temperature, °C 


x- and K, 
i i 


h i 


U 


thickness and thermal conductivity of the 
1 i ' th layer of the structural material 
inside convective heat transfer coefficient, 
kJ/m^ - h - °C 

2 

overall heat transfer coefficient, kJ/m -h-°C. 


In the present study, the values of 1 s are taken from [3], 
The inside convective heat transfer coefficient is found to 
be function of temperature difference ( A t = T , - T. ). 

W ,X X 
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For walls; 

h ± =4.25 (A T) 1 / 4 , kJ/h - m 2 - °C 2.121 

For ceiling: 

h ± = 3.48 (At) 1//4 , kJ/h - m 2 - °C 2.12(1 

The Newton-Raphson iterative procedure has been used to 
get the values of 1 h^ ' . Equation 2.9 is used to calculate 
the heat transfer through the structure; after considering 
time lag and decrement factors [2 ] . The effect of these 
factors on the cooling loads has been discussed, in 
Chapter 4 . 

2.7 HEAT TRANSMISSION THROUGH GLASS WINDOWS: 

The incident radiation may be reflected/ partly 
absorbed and the remainder transmitted through the glass 
materials. Thus, in general, it may be written as : 

r g + a g + ~ 1 (2.13 

where r , a and t are the reflectivity, absorptivity 
g g r 

and transmittivity , respectively. The equations for 
calculating these thermal properties for the glass windows 
are given in Appendix -D. 

From the basic concept of heat transfer, the rate 
of heat gain through the glass is given by. 
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g glass JJ s* t rD‘ 


Ip.+t .1 ,+t _.i_+h. (T . 
D rd d rR R x,g g,z 


- T. ) 
x 


( 2.14) 


where F g is the sunlit fraction of the window surface, 

^'rD' t rd / and t rR are the transmittiv * ties of the direct, 

diffused and reflected radiations, respectively . 1^,1^, 

and J R are the intensities of direct, diffused and 

reflected radiation. T . and T. are the temperatures of 

g,z x 

inside surface of the glass and indoor conditions, 
respectively. 

In reality, a glass window is not fully exposed 
to direct sunlight, but partially shaded. Referring to 
Fig. 2.2(a), it can be written as ; 


X = b . tan a 2.15(a) 

Y = b. tan 6 2.15(b) 

where the profile angle 6 is related to the sun's altitude 
angle p and wall solar azimuth angle a by. 


tan 6 p 


tan P 
Cos “ 


( 2.16) 


Finally, the value of sunlit fractions » is given by, 

s 

F = 1-P. .tan 6 - P„,tan a + P-,. P 0 . tan 6 tan a (2.17) 

s 1 P 2 12 p 

where , 

P^ = b/a and P 2 = b/c. 


m 
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Referring to Fig. 2. 2(b), the energy balance for the glass 
sheet can be •written as. 


P s* I D* a 


gD +I d* a gd + I R‘ a gR 


h. 

1 


( T . — T . ) + 

,g g.i i 


h (T -T 
o,g g,o o 


) 


( 2.18) 


where h- and h are the surface coefficients given by 
i,g o,g 

h Q ^ g = 3.96xlO" 3 + 9.89 2xl0" 4 .V Qy kJ/m 2 -°C -S 2.19(a) 

h, = 0 . 237 +( V . /Height ) ° * 3 , kJ/m 2 -°C-W 2.19(b) 

i /g 

where V Q is the outside air velocity in Km/h and is the 
indoor air velocity in m/h. 

Combining the equations (2.14) and (2.18),, the heat transfer 
through windows and doors including solar heat gain is 
calculated as 


°g = ^ F s * t rD* I D + fc ra • J d + + 


(F .a +a j.I.+ a I„) 

S QD D gd d CfR R - . rr / rn rn ) 1 

■ d + +u( ° 11 ]k 


( 2 . 20 ) 


where. 


u = 




h i,g 


o,g 


N = number of windows and doors. 
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The details of windows and doors are given in Appendix -E* 

The total heat transfer through the structure is given 
by the sum of equations (2.9) and (2.20) as: 

Q gt = 0 t + Q g (2.21) 


2.8 INFILTRATION LOAD; 

Infiltration is the leakage of outdoor air into 
a building through cracks and opening caused by pressure 
difference across the boundary surfaces. The net exchange 
of air may lead to both heat and moisture gain for the 
space. The volume of infiltration is related in terms of 
the room volume and taoular values are available in [3]. 
The expression for the infiltration load is; 


0. 

infil 


V 


room 


air 


x N. ( h 
ACri a 


h. ) /24,1c w 

i v a 


( 2 . 22 ) 


where _ is the volume of the room/ v . , is the 

room ' air 

specific volume of air at outside temperature, I<r ACH is 

the number of air changes/day, h and h. , are the enthalpies 

a -X 

of outside and inside air, respectively. 


2.9 VENTILATION LOAD; 

The introduction of outside air into the conditioned 
space is an important factor,, which is necessary to maintain 
oxygen and odour levels as per standard practice. Ventilation 
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codes and standards recommend a minimum rate of 30 cfrn/ 
person. But the recent analysis shows that it is more 
efficient at 5 cfm/person [19] on the basis of economic 
consideration. Since the present problem deals with 
economic aspect, this value has been adopted for the 
ventilation load calculation, which is given by 

• N 

“ventil “ - °°v Upan S * C Vent i 1/person) * ( h a - h ± > 
air 

( 2.23) 


where v .is the specific volume of air and h , h. „ are 
air a i,a 

the enthalpies of outside and inside air. 

2.10 MISCELLANEOUS LOAD: 

Apart from transmission, there is a considerable 
amount of heat energy generated inside the conditioned 
space from the occupants, lights, electrical equipments 
and appliances. The thermal energy release frcm human 
body at different degree of activities is taken from [l^i^L 
In the present analysis, the degree of activity is taken 
as the average cf all activities and the corresponding 
heat release has been calculated. The instant rate of 
heat gain from lights and other electrical appliances are 


calculated as. 
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Q , = (Humber of electrical appliances x rating of 

eiec the appliances) + (Number of lights x ratings) x 

cooling load factor 


(2.24) 


In the present work/ the cooling load factor is taken as 
0.88 [15] . So the miscellaneous load is calculated as 

Q = Q + Q , (2. 25) 

mi sc elec occupants 


The total cooling load is calculated from: 


Q total Q st + °infil + Q Ventil + Q misc 


( 2.26) 


Using this procedure/ in the present analysis, hourly 
cooling load calculations have been carried out for 
a building model. Fig. 2.3, for Kanpur. For this estimation, 
a generalized computer programme has been developed with 
great care and is given in Appendix-G. 
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CHAPTER 3 

PROBLEM FORMULATION AND OPTIMIZATION 

3 . 1 INTRODUCTION : 

The problem formulation for optimization is done 
on the ]»asis of total cooling load based on the design 
condition for the equipment selection. The total cost 
calculated on the basis of fixed and running costs for 
comfort airconditioning constitutes the objective function 
in terms of the operating variables having a set of 
constraints. The optimization problem can be stated as 

which maximizes or minimizes f(X) subject to 

-* 

the constraints g m (X)is O, m = 1/2/... where X is an 

n -dimensional vector called as the design vector/ f(X) 

is the objective function and g m (X) are the constraints . 

3.2 DESIGN VARIABLES: 

Any engineering system is described by a set of 
quantities some of which are viewed as variables during 
the design pricess. In general/ certain quantities are 
usually fixed at the outset/ known as preassigned parameters 
and all the other quantities are treated as variables in 
the design process and collectively represented by the 
vector X. Even though there are many variables which control 
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COMPRESSOR 



(b) 


Fig. 3.1 (a) Single-stage vapour-compression system 
(b) P-h diagram for single-stage vapour - 
compression cycle. 




Cl 


the airconditioning performance, inclusion of all these 
variables will lead to enormous computational effort. 
Here, an optimal prediction effort has been made to 
select the variables judiciously. 


3.3 SYSTEM DESCRIPTION: 

In the present study, a single-stage vapour 
compression refrigeration system is taken for the 
analysis having refrigerants R— 12 and R-22. The cycle 
is shown schematically in Fig.3<l* To simulate the cycle 
close to represent actual processes, a pressure drop in 
the compressor valvSs, heat rejection in the compression 
discharge line and the degrees of subcooling of condensate 
and superheating of evaporator vapour are considered. 

The pressure drops in the condense r, evaporator and 
piping are ignored. In the present analysis, air-cooled 
condenser has been incorporated in the system. 


In the design of the cycle, the pressure drop at 
the suction line is assumed as 0*2 toax [17 ]* The 
increase in pressure at the delivery of the compressor 
is taken as 0.4 bar [ 17] . In the present optimization 

analysis, four variables are considered. They are , 


X(l) 
X( 2) 


Degrees of rise in the condensing temperature, 

v° c 

Degrees of fall in evaporator temperature. 
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X( 3) : Degree of superheating/ T , °C 

s 

X( 4 ) : Degree of subcooling/ T c ,°C. 


3.4 ESTIMATION 01 OBJECTIVE FUNCTION: 

In the present study, the objective function is 
taken as the total cost for the comfort airconditioning 
system, which combines the initial and running costs of 
the system. If C^. and are the initial and running 
costs per year respectively, then the total cost/year 
( Cp) is given by 


where 


and 


C T 


C I “ 


comp 


evap 

C cond 

Si ower 
CL 


C^. 4 * Rs/ year 

C 4* C 4* C 4“ P 

comp evap cond blower 

= Cost of the compressor, Rs/year 

= Cost of the evaporator, Rs/year 

= Cost of the condenser, Rs/year 

= Cost of the blower , Rs/year 

= c 1 + c 2 + c 3 + c 4 

= Cost of the compressor power, Rs/year 

= Cost of the blower power, Rs/year 

= Cost of the fan power, Rs/year 

= Maintenance cost, Rs/year 


(3.1) 

3.1( a) 


3.1(b) 




The 'Present Worth Method' is used for the analysis of 


the economic model 
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3.5 OPTIMIZATION TECHNIQUE: 

The presence of constraints in a nonlinear progra- 
mming problem creates more problems while finding the 
minimum. Several situations can be identified depending 
on the effect of constraints on the objective function. 

The simplest situation is when the constraints do not 
have any influence on the minimum point. From the 
application point of view on the real problems, there is 
a possibility of having two or more local minima for 
the minimisation of the problem. This situation is entirely 
due to the nature of the objective function contours. By 
taking into account of all these possibilities, in the 
present analysis, one of the indirect methods of non— linear 
programming constrained optimization technique, known as 
'Interior Penalty Function 1 [ 16 J is used, 

3.5.1 THE PENALTY FUNCTION APPROACH : 

Penalty Function Methods transform the basic 
optimization problem into alternative formulations such 
that numerical solutions are sought by solving a sequence 
of unconstrained minimization problem* For analysis 
purposes, a basic optimization problem is taken such that 
X minimizes the objective function, f(X) subject to the 
constraints < 0, m = l,2,»..*nj* This problem 

is converted into a unconstrained minimization problem by: 
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0 k = 0(X, r k ) = f(X) + r k S 1 G. [ g(X) ] (3,2) 

j=l J 

-> 

where Gj is the function the constraints, g( X) , and r^ 
is the positive constant known as the ‘Penalty Parameter 1 , 
The function, 0(X, r^) , is optimised for a decreasing 
sequence values of the penalty parameter, r^. Evidently, 
if the 1 r^’ is veiry small, the optimization of 0(X, r^) 
amounts to optimization of the objective function, f(X). 

3,5.2 INTERIOR PENALTY FUNCTION METHOD: 

"*1 

In this method, the function Gj[g(x)J of the 
equation (3.2), is substituted by - — . The 

9m (X) 

unconstrained minima of 0(X, r^) , all lie in the feasible 

region and converge to the solution of f(X) as is 

varied in the decreasing sequence. So the function 
■+ 

0(X, r^) is defined as, 

n l 

0C X, r. ) = f(x) - r. S i — ( 3.3) 

k * m=l g m CX) 

Here, the values of function 0 will always be greater 

than f, since g ffi (X) is negative for all feasible points, 

-* . 

X. Since the unconstrained minimization of 0(X, r^) is 
to be achieved by changing the sequence of penalty 
parameter, r^, in the decreasing order, its initial value • 
should be chosen such that quick convergence is obtained. 
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Thus, for any feasible starting point, X^ , the value of 
r^ can be taken as. 


0.1 f (x x ) 


n l 

XT 

m=l 


r l S 


f (x x ) 


9m <V 


m=l 


C 3.4) 


9m tX l> 


Once, the initial value of r^. is chosen, the subsequent 
values of have to chosen such that 

2" s r* 3 • 5( a ) 

r k+l < k 

In this regard, for convenience purposes, it is taken as, 

r, - = D r, , where 0<D <1 3.5(b) 

k+1 k 

where in the present analysis, D has been taken to be 0.1* 

In the present analysis, the unconstrained 
minimi zation of the penalty function 0(X, r^) is done by 
'Davidon-Fletcher-Powells Method 1 [ 16] • This is the best 
general purpose optimization technique making use of 
the derivatives that is currently available. The algorithm 
for this iterative technique is as follows: 

(a) Start with an initial feasible point X^ and 
a nxn positive definite symmetric matrix, [HJ 
is taken as [i] • 


which 
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( b) Compute the gradient of the function, 
at the point and set s = - [Hi] . . Normalising 


this, we obtain, S 


the number of variables. 


.S 1 [(V0 1 ) 2 ] 1/2 

1=1 x 


, where 'n' is 


In the present formulation, central difference 
scheme has been adopted for gradient evaluation. 


(c) Find the optimal step length in the direction 


of S Ni and set X^ + ^ = (X i + T j_ s Ni ) • In the present 

★ 

study, the determination of T. is done by 'Cubic Inter- 


polation Method 1 [16]. 


< e , where e ^ is very 


(d) Test the new point X ±+1 for optimality using 
( V0 i + 1 “ V0 i ) 

the condition £ e 1 ' where e ^ is very 

small. In the present case, the termination criteria 
is e . < icf 3 . If this criteria is not satisfied [H] 

is updated as follows: 


[ H ±+1 ] 


[Hj+P ,] + P <1 


where , 


* T 

T . . S _ . . S*-, 

l Nl Ni 

<4 • °i> 


-(HA) 

<4V°i> 


Q ± = v 0 (X ±+1 ) ~v 0(x ± ) 
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START WITH AH INTERIOR FEASIBLE POINT lt{ 

~ T 



CALCULATE THE VALUES OF OBJ- 
ECTIVE FUNCTION AND CONS - 
[TRAINTS GG (M) 

I 4 

CALCULATE 0(X,R K ) W) -R K I^T/GG(M) 
HERE, R 1 IS TAKEN AS * CX OBJ* FUN /GG(M ) 
WHERE, OCCO 

— ±? 


FIND 0(VXp « V*()Tj ). SET Si * -[Hi] x V0(tfj) 
WHEN l«1, [H]«[I] 


|5ET )Tj+l »5fi »(Tj )sf. WHEN iO, T ■ 1 
| FI N D ftdfivjT 



YES 


Jnq 

| Ql * V0i+rVI»i 

TZ Z E ZZI 

FIND Tf BY CUBIC INTERPOLATION METHOD 

" -r 

|CALCULATE [Mj]* Tj* (Sj)(Sj r )/<Sj T ) (CM) | 

[Ni]*-(HfQi)(HiQi)VHi Q| t 

X 


[SET [Hi+i] = fHi]*[Mj]*[Ni] 
^ 


JSET 7j*=7 i+1 
tHj] -[Hj*,] 



SET Rk+I'^^K^ 


^OPTIMUM * ^i*1 

WHERE OCCO.O 


* optimum ® *1 


-~ t SET R K *= ft K*1 1 


Fig. 3.2 Flow chart for the interior penalty function method 
using Davidon - Fletcher - Powell's unconstrained 
optimisation technique. 





Fig. 3.3 Flow chart for cubic interpolation. 



35 


(e) Set the iteration number as i = i+l and 
repeat the procedure from the step 'b', until convergence 
is achieved. Flow charts for 'Interior Penalty Function 
Method' and 'Cubic Interpolation' are schematically presentf” 
in Figs. 3.2 and 3.3, respectively. 

3.6 ESTIMATION OF LOADS ON EACH COMPONENT OF THE SYSTEM: 

The components of a refrigeration system have to 
oe selected on the basis of the design conditions. The 
power requirements of the apparatus are evaluated for 
the design conditions and a factor of safety is multiplied 
to ensure the capacities of various components from 
the safety point of view. The condensing and evaporator 
temperatures of the vapour— compression systems are considered 
to be the sole factors which influence effectively the 
design of the airconditioning system. The evaporator 
temperature is selected on the basis of the indoor conditions 
whereas the condensing temperature is selected on the 
basis of design temperature. 

So the condenser temperature is taken as, 

T. = T + AT + X(l) (3.6) 

h o a 

the design outside air temperature/ °C 

the approach temperature of the condenser , °C 

degrees of rise in condensing temperature, 
taken as a variable, °C. 


wnere , 


A T 


a 

X( 1) 
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The evaporator temperature is chosen on the basis of 
inside room temperature. So it is taken as, 

T. = T. - AT - X(2) ( 3.7) 

1 1 a 

whe re, 

T^ = the inside room temperature, °C 
At, = the approach temperature, °C 

a 

X( 2) = the variation in tb\e lower temperature 
which is taken as a variable, °C. 

In the present analysis, the values of the approach 
temperature have beenpas 6°C, 10°C and 14 C. To justify 
these values, experiment has been done on an air conditioner . 
The pressure gauges were fitted in the suction and discharge 
line. Fig. 3.4. The airconditioner has been operated and 
temperatures of air leaving the condenser and evaporator 
were taken. The air was supplied to both at about 30°C. 

The average approach temperatures for the condenser and 
evaporator sides were found to be about 10°C. 

The value of T^ is also checked up in the basis 

of the apparatus dew point temperature, T adp . Because 

the lower temperature, T^, should alvays be less than 

T , . The T , is calculated in terrfls of effective sensible 
adp adp 

heat factor, ESHF , which is given by 
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ESHF 


1.026 (T. 

l 


T , ) 
— §^E_ 


1 .026 (T, ~ T , ) + 2410 (W. - W , ) 
1 adp i adp 


(3.8) 


where , 

ESHP = 

( ERSH + ERIri) 


ERSH = RSH + (BPF) . OASH 

ERLH = RLH + (BPS') . OALH 

= Specific humidity inside the room, kg/kg of dry a 

W adp ~ Specific humidity at apparatus dew point 
temperature, kg/kg of dry air. 


RSrf and RLH represent the room sensible and latent heats 
respectively. OASH and OALH denote the outside air sensible 
and latent heats, respectively. The effect of the bypass 
factor ( BPF ) on the performance of the cooling apparatus 
is discussed in [4"] and 0.15 is the value recommended for 
row coil. By using Newton-Raphson iterative procedure, 
the value of , is obtained, for the known values of 

ESHF. 


Referring to the Fig. 3.1, after taking the degree 
of superheat X(3), the temperature at the state point * l ‘ 
is gi ven by , 


T 1 = T x + X(3) - A T i ( 3.9) 

where At, is the decrease in the lower temperature before 
compression, °C. 
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Similarly, T 2 = T h + AT h 


(3.10) 


where A T^ is the increase in the condensing temperature 
after compression, °C. 

The enthalpy values, h pf the refrigerants, R-12 and R-22, 
have been computed from tne functional relations of the 
temperature of the state. These relations are given in 
the Appendix-F, 

Ftom the Fig. 3 .1 , = T^, - X( 4) (3 


where X(4) is the degree of subcooling, °C. 

* 

The mass flow rate of the refrigerants, m cef ' is found 


from : 


m 


ref 


Q >; SF 

c 

< h l - h 4> 


where Q q is the cooling load as given in Sec .2.10 and 
SF is the factor of safety taken as 1.25 [4]. 

The compressor power is given by: 


' comp 


m ref (h 2 - V 
T) 


, kw 


( 3.12) 


( 3.13) 


where r. is the mechanical efficiency of the compressor 
taken as 0.85 [12]. 

For the air-cooled condenser, the load is 


having 


°cond " “ref (h 2 ' h 3> ' kW 
h 2 = h l + ^ A h 'isentropit/ 11 o 


( 3.14) 
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The equation (3.13) can also be equated to the heat gain 
by the cooling air as 


Q 


cond 


- m . x c ( T - T . . ) , kW 

air p . o y c air,i 

o-L J- 


3 . 15( a) 


where , 


m 


air 


P • 
air 


T 


OyC 

T . 
air yi 


mass flow rate of cooling air, kg/s 
specific heat of air, kJ/kgi°C 

the outlet temperature of the air/°C 
the inlet temperature of the air,°C. 


In the above equation, the temperature differential is also 
equal to the variation in the condensing temperature (T^) 
denoted by X( 1 ) . So the equation 3.15(a) can be rewritten 
as. 


Q 


cond 


m air * c o * Xd), kW 
v air 


3.15(c) 


By combining equations 3.14 and 3.15b, 


m (fu 

ref 2 


h 3 > 


m 


' air 


c p . * 

*air 


x( l) 


/ kg/s 


( 3.16) 


From equation (3.16), the volume of the air handled by 
the fan for the condenser is calculated as. 


Vol . 
air 


m . x 60 _ 

, m 3 /min 

air 


(3.17) 


3 

The fan power is correlated with its capacity (m /min) 
in the form of a quadratic equation. For the purpose, a 




F - FAN 

S - WATER SUMP 


ROOM 


A -OUTDOOR DESIGN CONDITION 
AB - ADIABATIC HUMIDIFICATION 
BC-HEAT GAfN IN THE ROOM 
C -INDOOR DESIGN CONDITION 


26.2°C 30°C 41.2 # C 
DB OB DB 


Fig. 3.5 (a) Evaporative cooling arrangement in roo 
(b) Evaporative cooling process. 


set of volume flow rate and power consumption readings 
have been taken for three blowers , experimentally. 
These readings are plugged in the quadratic equation 
and they have been solved for the constants. So this 
procedure finally results in the form: 


blower 


3.429 (Vol , )* + 1.143 (Vol . J + 37.67 

aJ- X d.JL L. 


t kW 


looo x r) 


m 


( 3.18) 


where , 


Vol 

r) 


air 

m 


3 

The volume flow rate of air, m /s 

The mechanical efficiency of the fan (0.85) 


3.7 EVAPORATIVE COOLING PROCESS: 

Figure 3.5(a) exhibits the evaporative cooling 
arrangement for the room. The unsaturated outside air 
comes in contact with the wetted surface where the evapo- 
ration of the water lowers the humidified air temperature. 
This process is schematically represented by process AB 
in 3.5 (b). The temperature of the humidified air is 
related to the desert cooler efficiency C 7 ^) as 

T dg = T db " r 'D ^ T db " T wb^ (3.19 

where, 

= dry bulb temperature of the outside air,°C 

T i = wet bulb temperature of the outside air, °C 

T, s= dry bulb temperature of the air leaving 
g the desert cooler. °C . 


43 


The volume of humidified air for the airconditioning is 
calculated from: 


d,air 


Q c 

^air ^a,o h a,i^ 


0 . 


air * 1 ‘ 026 • (T i - T dg } 


( 3.20) 


where P a ^ r is the density of the humid air. 


3.8 CONSTRAINTS FOR THE PROB LEM : 

As mentioned earlier, four variables, being restricted 
between some practical values, are taken into consideration 
for the optimization. So, these constraints are mathematically 


^pressed 

as 

follow 

s: 



4.0 

< 

x( 1 ) 

< 

15.0 

3 *21( a) 

2.0 

< 

X( 2) 

< 

5.0 

3.21( b) 

1 .0 

< 

X( 3) 

< 

4.0 

3 . 21( c) 

1.0 

< 

X( 4 ) 

< 

4 .0 

3.21(d) 


The constraints are transformed into equality constraints 
by taking the lower and upper limits, separately. As 
a result, they are transformed in the following form: 

'i 

/if 


GG( 1) 

= 4.0-XCl) 

3.22(a) | 

GG( 2) 

= X(l) - 15.0 

3.22(b) 

GG( 3) 

= 2.0 - X( 2) 

3.22(c) § 

GG( 4 ) 

= X( 2) - 5.0 

3.22( d) M 
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GG( 5) 

= 1.0 - 

■ X( 3) 

3 . 22( e) 

GG( 6 ) 

= X( 3) 

- 4.0 

3.22(f) 

GG( 7 ) 

= 1.0 - 

• X( 4) 

3.22(g) 

GG( 8) 

= X( 4 ) 

- 4.0 

3.22(h) 


After normalising the equations, the constraints are 
considered for the optimization of the problem. 


3.9 ESTIMATION OF INITIAL COST: 

As mentioned in Sec. 3.6, the design cooling load 
is calculated for Kanpur, in May as in [9] . After evalu- 
ating these loads, the equipments are selected on the basis 
of their capacities. In this regard, a due consideration 
has been given to purchase the equipments on the basis of 
recommended capacity ranges available in the market. A 
stepwise variation is adopted for the costs and capacities 
of the various equipments. This feature is very well 
incorporated in the computer programme. 

Thus, the total initial cost for the system is 
expressed as 


Cl - ( C com p + c eva p + C cond + ^lower^ X 

( 3.23) 


where. 


C = cost of the compressor, Rs 

comp 

C = cost of the evaporator, Rs 

evap 
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C , = cost of the condenser/ Rs 

cond 

C blower = cost the blower, Rs 

MF = Multiplication factor to take care of 
the miscellaneous cost (1-.2) 

L> s = life of the system, years. 

of 

When a hybrid system ^air conditioning and desert cooler 
is adopted between the months April and June, if feasible, 
the initial cost, Cj, is considered to be higher and 
given by, 

c l = c l + C Desert //L d (3 * 

where, 

Cj = initial cost of the airconditioning system, 
1 Rs/year 

C = cost of the desert cooler, Rs 

Desert 

= life of the desert cooler , years. 

Jn the present work, both these cases are studied 
separately and results are compared. 



3.10 ESTIMATION OF RUNNING COST: 

For the running cost, the electricity cost is the 
primary factor to be considered. Keeping in view, the 
variation in electricity cost, a futuristic statistical 
approach has been adopted to make the whole formulation 
as a generalised algorithm [13 ] . Proceeding as in [13], 
the electric charge ,( C e ) ,is found to be: 


/ 


24) 
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1.293 

0.97986 + exp (-0.09338 I) 


( 3.25) 


where/ I refers to the duration of the years from the 
present one. 


The effective electricity cost is given by [l5]i 

, L c 


°E = i 2 
1=1 


( 1+R) 


i-1 


( 3.26) 


where / 

C e = cost of electricity in the ith year,R S /kWh 
R = interest rate 

L = life of the machineries, years 

So, the running cost for compressor (c^) is given by. 


C-, = P x C„, Rs/h 

1 comp E 


(3.27) 


For the evaporator fan, 


C 2 = han * V Rs/h 


For the condenser fan/ 


C 0 = P. x C„, Rs/h 
3 fan E 


( 3.28) 


( 3.29) 


The running cost of all these accessories is calculated 
on the hourly basis by taking into account by sinusoidal 
variation of outside air temperature and summed upto get 
the running cost/year. The maintenance cost is taken 
to be 10% of the initial investment [15]. 





Mathematical ly, it can be expressed as , 


CL = (4.1) C - -- --- — “=■ — 

4 1 R( 1+R) L - 1 


(3.30) 


The total running cost is given by: 


°R = C 1 + C 2 + C 3 + C 4 / Rs /^ ear 


( 3.31) 


In the case of combined operation of airconditioning and 


desert cooler the running cost includes for the desert 


cooler only for the feasible period during April to June, 


Then the total running cost becomes/ 


°R “ C E * D - C rating * + C 4 


( 3,32) 


where CL, the effective electricity cost, Rs/kwh 
E 


D.C . The desert cooler rating/ kW. 
rating 


Thus, the total cost of the system, Q-p,is given by 


C T = C R + C I/ Rs/year 


(3.33)?:;; 


*i n Eq.( 3.30) when desert cooler is also used, is replaced $ 
by Cj of Eg. ( 3.24) . | 


"Thin operative cost is for the period during which it 
operates alone. 
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In the case of hybrid system/ it is 


C T = C R + Cj. / Rs/year 


( 3.34) 


For generalisation purpose/ the equations 3.. 33 and 3.34 
are modified as. 



( C 


R 


+ c I ) 


E 


kWh/year 


3.35C a) 




C 


E 


( c R + c ;> 



kWh/year 


3.35(b) 


Thus the ratio Qp/C^ can be written in functional form as, 

Qp/Cg = f(Xp, X , X , X 4 ) , where X ± ' s are the variables 
taken for the present study to get the optimum values for 
the minimum of the ratio Cp/C^. 

In the present analysis, all these possible combinations 
are incorporated in a generalised computer programme for 
optimization and the same is incorporated in Appendix— G. 
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CHAPTER 4 


RESULTS AMD DISCUSSION 


This chapter deals with the calculation of cooling 
loads on the basis of the theoretical and actual temperature 
variations having the effect of heat capacity and time lag 
factor of the structures. The economic criterion has been 
adopted for the optimum choice of the airconditioning systems 
using R-12 and R-22 refrigerants. A comparison of a hybrid 
system (evaporative cooling and airconditioning systems) 
with a conventional airconditioning system is discussed 
under optimum conditions. For estimation of the optimum 
cost/ Qp/Cg, the life of the airconditioner and the desert 
cooler are taken as 20 years and 10 years, respectively* 

The following are the results obtained in the present work. 



4.1 THEORETICAL AND ACTUAL TEMPERATURE VARIATION; 

Figure 4.1 gives the actual temperature variation 
with time as well as the suggested sinusoidal variation in 
temperature based on maximum and minimum temperatures of 
the day. The actual temperature is slightly higher from 
8 A.M. to 3 P.fl. than the theoretical value before the peak 
hour. Thereafter, the actual temperature is lower than 
that of the theoretical value. The areas under both curves 
reveal small difference ( about 4 to 7%) . For the airconditioning 
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system, the working hours have been taken to be 8 A.M. to 
* 

5 P.M. . An underestimation of cooling loads by using 
sinusoidal variation is found in this. This deviation is 
calculated for three years, 1982, 1983, 1985 and the same 
are tabulated in Tables 4*la to d, for different indoor 
conditions. From these values, the average deviation is 
calculated as 6%. Sc a multiplier 1.06 can be used to 
get the actual total running loads for 10 hours. 

The variation in the sol-air temperatures for 
different directions is shown in Fig. 4.2. It indicates 
that the maximum sol-air temperature occurs on the western 
facing wall at 3 P.M. 

4.2 EFFECT OF HEAT CAPACITY AND TIMELAG FACTOR ON THE 
COOLING LOADS: 

Figures 4.3(a) and 4.3(b) represent the effect of 
the heat capacity and timelag of the structures on the 
hourly cooling load estimation. From this, it is found out 
that the estimation of the cooling load without taking into 
account of heat capacity of the structures leads to an 
overestimation of 20.29% than that of the cooling load 


The computer programme is quite general and can be used 
for any duration. 



Time , hours 

Fig, 4.1 Variation in outside air temperature with respect to time in May at Kanpur. 






Time, hours 

Fig. 4.2 Variation in sol -air temperature and outside air temperature w.r.t. time at Kanpur. 
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calculated for the structure with the heat capacity. This 
difference is quite pronounce. When the wall thickness 
was halved , the above is found to be 18.02%. 
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4.3 OPTIMUM RESULTS FOR VARIOUS INDOOR CONDITIONS: 

For the indoor conditions of 24 °C DB and 65% RH , 
the design cooling load is 9,8% and 24.5% higher than 
that for the indoor conditions of 27. 5° C, DB , 56% RH 
and 30°C DB , 60% RH, respectively. Figures 4.4 and 4.5 
represent the variations in initial, running and total 
costs with respect to the rise in condensing temperature 
(T^/ of the system, for R-22 and R-12 refrigerants, 
respectively* In this case, the approach temperature is 
taken as 10° C, while the other design variables, the fall 
in lower temperature AT^ the degree of superheat A T gup , 
and the degree of subcooling AT gub are maintained at 
the optimum values. The optimum total cost / Qp/Cg, for 
R-12 refrigerant system is 4.752% less than that of R-22 
rigerant system for the same operating variables. 

Figure 4.6 shows that the optimum total cost, 

( Qp/c ), for a hybrid system of evaporative cooling and R-12 
air conditioning system yields 4% gain over the conventional 
airconditioner used from April to October, every year. This 
is applicable for the indoor conduction of 27.5°c dry-bulb 
temperature and 56% RH, with approach temperature as 10° C. 
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Figure 4,7 depiots the comparison of the optimum 
total cost between a conventional airconditioning system 
and the system with evaporative cooling. The optimum total 
cost for the conventional airconditioning system is 5.562% 
higher than that the system used with evaporative cooling 
facility. The variation in the optimum total cost (Gp/Gg) 
with respect to the fall in the lower temperature of the 
system, At^ is displayed in Fig. 4. 8, The effect of AT^ 
on the optimum total cost is less than that of A T^. 

Figures 4.9(a) and (b) show that the degrees of superheat 
and subcooling have relatively little effect on the optimum 
total cost of the system. 

Figure 4.10 represents the variation in the 
optimum total cost with respect to the approach temperature 
of the systems for various indoor conditions. The approach, 
temperature is almost proportional to the optimum total 
cost or vice versa. The various values of optimum variables 
for R-12 and R-22 refrigerant systems are given in the 
Table 4,2, along with the optimum total cost, penalty function 
values and the penalty parameters. 
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TABLE : 4.2 : (CONTINUED) 



t/i C*/Ce, kWh/Yr, 





Fig. 4.4 Variation in cost with rise in condensing temperature 










10 3 l x10 3 l xIO 3 


For indoor conditions of 27.5°C DB and 56%RH 

Refrigerant R-12 
^approach 5 1Q°C 

SYSTEM WITHOUT DESERTCOOLER 

SYSTEM WITH DESERTCOOLER 


7.5L 3.8L 3.2L 


0 4 8 12 16 
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Fig. 4.6 Variation in cost with rise in condensing temperature 
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Fig. 4.8 Variation in cost with fall in evaporating temperature. 
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CHAPTER 5 

CONCLUSIONS AMD SUGGEST IONS 


5.1 CONCLUSIONS: 

i) A multiplication factor has been found out for 
the estimation of actual cooling load with 
the help of maximum and minimum temperatures 
of a given location. 

ii) The heat capacity and timelag factor render the 
cooling load about 20% less than that of 

without these factors. 

iii) The evaporative cooling system can be substituted 

for the conventional air-conditioner during the 

period of April to June in the case of inside 

design condition of T^ b = 30°C/ 0 = 60% and 

T = 27 5°C/ 0 = 56%. 
do 

iv) The implementation of the hybrid evaporative 
cooling and airconditioning system is found to 
be more economical than the conventional 
airconditioner , under optimum conditions. 

v) For the inside conditions of 27. 5° C D3, 56% RH 
and 3 0°C DB, 60% RH, the suggested hybrid 
system yields a gain of 4 % and 5.6% over the 
conventional airconditioning system, 
respectively • 
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vi) The optimum total cost of the comfort 

airconditioning using R-12 is cheaper than 
that of R— 2 2 system by 4.7 5% for the inside 
design conditions of 24°C D3 and 65% RH . 

vii) The variations in condensing and evaporator 

temperatures have more influence on the total 
cost of the system than the degrees of sub- 
cooling and superheating. 

viii) The variation in total cost of the airconditioning 
system is almost proportional to the approach 
temperature. 

5*2 SUGGESTIONS: 

i) a combination of evaporative cooling and conven- 
tional airconditioner can be fabricated as a 
single unit. The overall cost can be estimated 
on the basis of actual cost for commercial 
venture. 

ii) The cost of the building may also be included 

in the optimization in order to make the approach 
even more general* 

iii) This approach can also be extended to the 
vapour-absorption system under optimum conditions* 


f 
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APPENDIX-3 
Table -3 

VALUES FOR SOLAR INTENSITY CALCULATIONS 


■ — ~ r m 

i 

Month • ! 

i 

i 

i 

i 

i 

Equation of 
t i me/mi nutes 

~i 

i 

| Declination 

J of angle 
! degrees 

i 

j A 1 

! W/m 2 

1 3 1 

1 

t 

t 

ml 

C 1 

Jan. 

-11.2 

-20.0 

1230 

0.142 

0.058 

Feb. 

-13.9 

-10.8 

1214 

0.144 

0.060 

March 

- 7.5 

0.0 

1185 

0.156 

0.071 

April 

1.1 

11.6 

1135 

0.180 

0.097 

May 

3.3 

20.0 

1103 

0.196 

0.121 

June 

-1.4 

23.45 

1088 

0.205 

0.134 

July 

—6 . 2 

20.6 

1085 

0.207 

0.136 

August 

-2.4 

12.3 

1107 

0.201 

0.122 

September 

7.5 

0.0 

1151 

0.177 

0.092 

October 

15.4 

-10.6 

1192 

0.160 

0.073 

November 

13.8 

-19.8 

1220 

0.149 

0.063 

December 

1.6 

-23.45 

1233 

0.142 

0.057 


A x = Apparant solar irradiation at air mass = 0 , W/m 2 
B 1 = Atmospheric extinction coefficient 
C = Diffuse radiation factor. 
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APPENDIX - C 

EQUATIO NS FOR TH5 3STIMATI0N OF PRIMARY ANGL5S 

The equation for^LST, local solar time 

LST = LCT + Equation of time (C-l) 

where^ 

LCT = 1ST + (82.5-Longitude of the 

location )/l 5 y hrs (c-2) 

The equations for the hour angle are given by 

a h = (LST - 12.0) .15 / for LST < 12.0 

(C-3) 

= (12.0 - LST) .15, for LST > 12.0 

The relationships between the latitude angle 1, the hour 
angle declination angle d, altitude angle p , and 

azimuth angle y are , 

Sin P =( Cos 1. Cos h. Cos d + Sin 1. Sin d) (C-4) 

Cos y= Sin (Cos 1. Sind - Cos d. Sin 1. Cos a^) ( C-5) 

For wall solar azimuth angle a t 

cc = -y 1, for east and west facing walls, 

= y , for north facing walls, 

= I (Y _7t) I , for south facing walls, (C-6) 
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The relationship for the incidence angle , © is 

cose = (COSI3) ( Cosa) (cosi^ + ( Sin ) (sins*) for 
tilted surfaces, with tilt angle jZS^ 

= (CosP) (Cosa), for vertical surfaces 
= Sirp,- for horizontal surfaces. 
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APPENDIX - D 
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APPEND XX -D 

EQUATIONS FOR THERMAL PROPERTIES OF GLASS 

The relationships for the thermal properties of the 
glass / t ran smith, vity^t^, reflectivity p. ^nd absorptivit y f a t 


t r = (1-r') 2 . a c /(l - r^. a 2 ) 


( D— 1 ) 


, 1 v 2 2 
( 1 - r* ) . a_ 


r = r‘ + 


9 (l-r' 2 .a 2 ) 


(D-2) 


a = l _ r ‘ 

g 


(1 - r )* . a. 


(1 - r* a c ) 


( D— 3) 


where , 


Sin 2 (0 -©') 


i. r a 

9 L 9 1 

Sin {Q + © ) 

g g 0 

Sin^ 6 

-E x /( 1 - 2^ 

P PT / v o / 

exp y n 


tan (© - 

g 

tan ( © •* 


where © Angle of incidence 

Of 

©g Angle of refracted rays 

E ■ Extinction coefficient 
c 

x Thickness of the glass (m) 


n Refraction index. 
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APP END IX-E 


building particulars for cooling load calculations 

Room Size = 8.4 x 4.0 x 3 .8, m 3 

Thickness of the wall = 0.4 m 

Thickness of the roof = 0.2 m 

Number of windows: 

On the Northern wall = 2 


On the Southern wall = 2 
On the Eastern wall = O 
On the Western wall = 0 

2 

Area of the windows = 1.08 m 

Number of doors: 

On Western walls = 1 
On the other sides = 0 
Number of occupants = 25 


Number of fans inside the room = 4 
Number of lights inside the room = 8 
Ratings: 

Fans = lOO watts 


Lights 


40 watts.. 
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APPUMDIX-F 


PROPERTIES OF REFRlGERAxrr.q 


Refrigerant R-12: 

c p( T ) = 0.59524 + 0 .001817 15 (T+15) 

S f ( T ) = 0.142093 + 0.33747 (T/100) - 0.0399606( T/100) 2 

+ 0.019868( T/100) 3 + 0.014023 (T/100) 4 
hg(T) = 188.86 + O,440278(T) - 7 ,02007( T/lOO) 2 

- 5.07651 (T/100) 3 _ 3.62545 (T/100) 4 

h f(T> = 36.1554 + 0.9 2bl08(T) + 6.89916 (T/100) 2 

+ 3.73414 (T/lOO) 3 + 5.91673 (T/lOO) 4 


Ref rig rant R- 2 2 : 

c p(T) = 0.60 + 0.00 0 5246 (T+4 5+0. 04 3686 (T+45) 2 )), 
for -60 < T < 10 

= 0.70114 + 0.0029529 ( T-IO+0. 01047 2( T-IO) 2 ) 
for 10 < T < 50 

S f (T) = 0.181332 + 0.43667 (T/lOO) - 0.0596756 (T/100) 2 

- 0.0234343 (T/100) 3 + 0.0638047 (T/100) 4 

h (T) = 251.06 + 35.4577 (T/lOO) - 19.4993 (T/100) 2 
9 

- 7.62005 (T/100) 3 - 11.6756 (T/lOO) 4 

h f (T) = 46.2102 + 1.20394 (T) + 6.828 (T/lOO) 2 

- 8.99387 (T/100) 3 + 20.0612 ( T/lOO) 4 

where y T in • °C, c p in kJ/k<g.°C, s f in kJ/kg. e K / h ,h f in kJ/kg . 
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